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Abstract 

The problem of motion of a single electron interacting with a periodic lattice of two- 
level systems is investigated within a spinless fermion model. The Green's function 
is calculated in a single-site dynamical coherent potential approximation which is 
equivalent to DMFT. The picture of one-electron density of states is obtained for 
various values of the tunnel splitting and coupling between the electron and two- 
level system. The occurrence of a band splitting with increase of the coupling is 
demonstrated. 
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The problem of interaction of current carriers with local excitations is one of 
the classical problem of solid state physics [1]. Recently, the interest in this 
topic has been revived in connection with investigation of highly correlated 
electron systems, in particular high-Tc superconductors where anharmonic lat- 
tice vibrations (i.e., those of apical oxygen) are assumed sometimes to play an 
important role [2,3,4,5,6]. The problem is relevant for both metals and semi- 
conductors where the electron spectrum is essentially modified by the influence 
of local excitations. 



Formally the local excitations can be described in terms of multi-level spec- 
trum of a strongly anharmonic system. In the simplest case one can use a 
model of two-level system (TLS) which is formally described by pseudospin 
formalism. This model is widely used for two-level systems in metallic glasses 
[7,8,9,10], crystalline-electric field (CEF) excitations [11], dissipative systems 
with tunneling states [12]. 
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In the present Letter we investigate the onc-clcctron spectrum in this model 
by using well established method — single-site dynamical coherent potential 
approximation (DCPA) [13,14,15,16] which becomes equivalent to dynamical 
mean-field theory (DMFT) [17] in the limit of infinite dimensionality of the 
lattice (d — > oo, which corresponds to Gaussian bare density of states for hy- 
percubic lattice) or infinite nearest neighbor number {z oo, Bcthc lattice). 
In those limits the equivalence of both methods was demonstrated explicitly 
in [18]. From the practical point of view, DCPA difi^ers from DMFT by that 
the former retains the shape of the bare density of states for finite d. Both 
the methods enable us to map a rather complicated TLS lattice problem onto 
a single-impurity problem which is readily solvable for one electron in empty 
conduction band. 

We consider the spinless fermion model describing fermion-pseudospin inter- 
action: 

H = Y1 tij^i^j - IZ(7iTi)^i - I^(hiTi) (1) 

ij i i 



where q, q are fermion creation and annihilation operators on a site i; tij 
are transfer integrals, rii = clci, Ti are pseudospin-1/2 operators; 7^ is the 
vector of coupling between fermions and pseudospins, hj is pseudomagnetic 
field vector, the tunneling frequency being = Aj = ^^(hjhi). 

Further on we consider translationally invariant case where site dependence 
of both coupling constants and tunnelling frequency is absent 

^ = E ekcLck - A E(t+ + r-)n, - /i. E (2) 

k i i 



with ek being electron band spectrum, and we have put 7^ = after pseudospin 
quantization axis rotation. 

Note that the model (1,2) can be mapped onto various models in different re- 
gions of its parameter space. For example, at 7^ = (0, 0, 7^) and hj = (0, 0, h^) 
we get the Falicov-Kimball model [19]. For large pseudospin values we ob- 
tain by using the Holstein-Primakoff representation the Holstein model [20,21] 
which is widely used to describe small polaron formation, in particular in 
molecular crystals. 

We calculate the one-particle Green's function in the case of single conduction 
electron, 

Gk(^) = ((Ck|4)). = [z-e^- Sk(^)]-^ (3) 
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where Tik{z) is the electron self-energy. 

In the local self-consistent approximation (which corresponds to dynamical 
CPA and DMFT) the quantity Sk(2;) is replaced by the momentum-independent 
local self-energy, Ek(2;) — > '^ioc{z) and the expression for the on-site Green's 
function G{z) — J2k Gk{z) takes the form 

G{z)^Ro{z-Ei,,{z)), R,{z)^^^—. (4) 

k ^ Ck 



The local self-energy is obtained from the solution of the auxiliary single- 
impurity problem, the corresponding model parameters being determined from 
the self-consistency condition G{z) — Gioc{z)- The local Green's function is 
given by 

- Sloe (2;) (5) 



where R^q'^{z) is the resolvent of the single-impurity problem, with the inter- 
action at the impurity site being switched off. We obtain the self-consistency 
condition as 



Introducing F{z) by F{Ro{z)) = z to exclude T,\oc{z) from the expression (6), 
one can write down the self-consistency condition in a more familiar form 

[R'^%z)]-' = z-F{G{z)) + G-\z), (7) 



if one prefers a DMFT-like formulation (see [17]). 

Further on we solve the auxiliary single-impurity problem with the Hamilto- 
nian on the fictitious one-dimensional lattice 



n=0 n=0 

-Xclco{r+ + r-)-h,r' (8) 

by using the method described in [22] to find Gioc{z) for a single current carrier. 
Here e„ and en,n+i are on-site energy levels and transfer integrals, respectively, 
for that fictitious lattice, n enumerates the lattice sites. Then we write down 
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the on-site impurity scattering matrix for the pseudospin projection a, 

t (z) - A^<-(^ - ^A) .9^ 



and introduce the average T-matrix 



t{z)=J2PJa{z): (10) 

P, = ^[l + atanh(/3A/2)], A = /i„ /3 = ^ 
to obtain the standard expression 

= 1 + Rrl)t{z) • ^^^^ 



As follows from the above treatment, Gioc(^) can be considered as a solution to 
the problem of scattering by the random substitutional impurity in the lattice 
with the bare resolvent Rioc{z). In such a case, the quantity A^-Rg°'^(z — aA) has 
the meaning of the random dynamical scattering potential for the impurity 
distribution P„ 



a- 



For the lattice model the quantity Sioc(^) calculated under the assumption 
R\oc{z) = Ro{z) corresponds to the average single-site T-matrix approxima- 
tion (ATA) , and the fully self-consistent solution corresponds to CPA for the 
disordered alloy problem [23]. This analogy enables us to use bonding and 
anti-bonding state classification. 

We use the semielliptic bare density of states (the Bethe lattice with infinite 
nearest-neighbor number which corresponds to DMFT situation), 

m = -^vd^, (12) 



D being bare half-bandwidth, and organize our numerical procedure as fol- 
lows. For a given R]^^{z) we calculate T,ioc{z) from (11), and then G{z) from 
(3,4). Using the self-consistency condition in the form (6) or (7) we recompute 
updated R^q'^{z) for the next step of the iteration procedure. The iteration 
process rapidly converges in a few steps. We also calculate Gioc{z) from (5), 
but for the first iteration only. 

As an initial condition for the iteration process we choose R^q'^Iz) = Rq{z). In 
that situation Giodz) obtained at the first iteration gives us an exact solution 
of the single- impurity problem on a lattice with the bare resolvent Ro{z) (i.e.. 
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Fig. 1. Density of states for the semielliptic bare band with D = 1, pscudomagnetic 
field A = 0.02 and different values of coupling parameter A in three approximations 
(a) single impurity in a lattice, left (b) dynamical ATA, center (c) dynamical CPA, 
right 

the electron-pseudospin interaction is switched off on all the sites but one). 

As for the first iteration lattice G{z), it has clear meaning of that calculated 
within the average single-site T-matrix approximation, speaking in an alloy 
analogy language. 

Basing on the above consideration, we discuss the following three approxi- 
mations (a) single impurity site in the lattice (b) dynamical ATA and (c) 
dynamical CPA. 

The results of calculations for zero temperature are shown in Figs. 1-2. One can 
see that with increasing A formation of bonding and anti-bonding states takes 
place, which is clearly marked in the cases (a) and (b). It should be noted 
that, unlike the scattering by static impurity potential (where bonding and 
anti-bonding states occur depending on the potential sign) , dynamical nature 
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Fig. 2. The same data as in Fig.l for A = 0.25 
of scattering in our model results in coexistence of both the state types. 

In the case (a) and for < A < D, the critical values Ac,b and Ac,a of the 
coupling parameter A for the occurrence of bonding and anti-bonding states 
are given by 



Ac,a=[i?o(^t)i?o(^t + A)]-^/^ 



Ac,b — 



(13) 



Here E-^ and Et are lower and upper bare conduction band edges, respectively, 
and we have taken into account that the quantities Imi?o(-Sb), ^^Ro{Eh ~ A), 
lmi?o(-E't) and lmi?o(-E't + A) are equal to zero. Provided that the bare DOS 
has a property 



N{Ei, + D + e)^ N{E, -D-e), 



(14) 
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(e.g, bare DOS is symmetric) the following equalities hold 



Rei?o(£^b- A): 
Rei?o(^b) 



-ReRo{Et), 



we have 



Ac,a — -^Cjb) 



(15) 



(16) 



and hence in this case bonding and anti-bonding states appear in the spectrum 
simultaneously. For semielliptical bare DOS (12) we put 



(17) 



to obtain 



A. 



Ac,b — TT 



1/2 



(18) 



In the dynamical ATA (case (b)) and for < A < we derive 
Ac.b = W — „ . „ s , An 



/2L) + A 



Rq{E^) ' 



Ro{Et 



(19) 



where the conditions Imi?o(-Eb) = 0, lmRo{Et) = 0, Rq{Ei,) < and i?o(-^t) > 
have been used. Here Ac,b and Ac,a are critical couplings for the situa- 
tion where bonding and anti-bonding state subbands have been already fully 
formed (e.g., the subbands are decoupled from a resonant state subband). 
Provided that the condition (15) takes place, one has 



Ac,b ^ '^c 



(20) 



and a bonding state subband appears in the spectrum at lower value of A being 
compared with that for anti-bonding state subband formation. For semiellip- 
tical bare DOS (12,17) we obtain 

Ac,b = i^^l-^, A,,, = i?^l + A. (21) 

For small A (Fig.l) the case (c) differs from the case (b) by a more smooth 
DOS behavior, which is a result of the full self-consistency in this approach. For 
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large A we obtain in DCPA a quite rich structure owing to intrinsic impurity 
dynamics at A 7^ 0, along with the standard band splitting. One can see from 
Figs. 1,2 that DCPA describes the formation of the energy gap in the band 
centre rather adequately, whereas dynamical ATA yields resonance states in 
the "pseudogap" even at large A. Prom analogy with disordered binary alloys 
one can assume that such a behavior is owing to lack of self-consistency in 
dynamical ATA (cf. the consideration of Hubbard-III-like approximation in 
strongly correlated systems [24]). At the same time, dynamical ATA gives a 
reasonable picture of bonding and anti-bonding states formation. 

The non-zero values of the critical values (13,19) are connected with the 
square-root energy behavior of the bare DOS at the band edges. One can 
expect that in the case of a rectangle- like bare DOS (e.g., two-dimensional 
case which may be also qualitatively described by DMFT) the bonding and 
anti-bonding states will occur at arbitrarily small A. A similar situation near 
the Permi level is expected in the case of finite band filling in the rigid-band 
approximation. 

The effective mass renormalization in our model is due to band narrowing. We 
obtain 



m* 9ReEi°^(E) 



m dE 



< 2. (22) 

E=E* 



Thus in our two-level situation the bandwidth renormalization does not exceed 
the value of two, and at large A only two subbands survive in the spectrum. At 
the same time, in the standard problem of small phonon polaron (which was 
also considered in DMPT [25] ) increasing electron-phonon interaction results 
in many-fold splitting of electron band. 

To conclude, numerical results even in the simplest particular case of our model 
demonstrate a rather rich and non-trivial spectrum picture. Investigations of 
the general model (1,2) seem to be of interest. 

The work was supported in part by Grant No. 4640.2006.2 (Support of Scien- 
tific School) from the Russian Basic Research Foundation. 
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